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We discuss the role of rare fluctuation effects in quantum condensed matter systems. 
In particular, we present recent numerical results of the effect of resonant states in 
Anderson's original model of electron localization. We find that such resonances give rise 
to anomalous behavior of cigenstates not just far in the Lifshitz tail, but rather for a 
substantial fraction of eigenstates, especially for intermediate disorder. The anomalous 
behavior includes non-analyticity in various properties as a characteristic. The effect of 
dimensionality on the singularity, which is present in all dimensions, is described, and 
the behavior for bounded and unbounded disorder is contrasted. 
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1. Introduction: A Brief History of Large Disorder and Rare 
Fluctuations 

Rare fluctuation phenomena received little attention in the early days of Solid State 
Physics. Ignoring them seems to have been the norm - what effects could such 
low probability phenomena possibly have, especially for a thermodynamically large 
system? On the contrary, the phenomenal success of Bloch's theorem and its ap- 
plication to wave phenomena in solids (e.g. electron energy bands, phonon disper- 
sions), led to a kind of "Blochitis" among the practitioners of the field. Materials 
not having proper translational symmetry (i.e. disordered materials like amorphous 
solids, glasses, liquids) were addressed using perturbative techniques starting from 
the uniform system, usually a crystal with an appropriate coordination number, 
and then averaging the variables representing the disorder. This led to the growth 
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of various mean-field like schemes for treating disorder, like the Coherent Poten- 
tial Approximation, the averaged T-matrix approximation, the Effective Medium 
Approximation etc. 

However, none of these methods quite captured effects that were primarily 
present because of disorder, as exemplified in the seminal paper by AndersorP in 
1958 entitled "Absence of Diffusion in Certain Random Lattices" . Anderson showed, 
for the first time, the existence of localized states in a disordered one-electron sys- 
tem, localized around disorder-specific locations in the system, in contrast to the 
Bloch waves in crystalline materials, which extended in a periodic fashion over the 
whole macroscopic system at all allowed energies. In our three dimensional world, 
localized states existed for moderate to large disorder. Further, as emphasized by 
MottPH these localized states were separated in energy from extended states; the 
separatrix was called the mobility edge; in recent years, it is also referred to as the 
critical energy of the localization-delocalization transition, Eq- Thus, with increas- 
ing disorder, localized states appeared first in regions of low density of states (band 
tails), and gradually the mobility edge(s) moved towards the center of the band, 
until at a critical disorder, all states became localized and no extended Bloch-like 
states were left. (In reality, extended states for moderate disorder also exhibit large- 
scale fluctuations that are disorder specific, and arc known in the field as mesoscopic 
fluctuations) . 

Soon after Anderson's paper, Mott and Twose^ showed that all electronic states 
were localized in one dimension for arbitrarily weak disorder (and hence Bloch's 
ideas broke down immediately upon introduction of disorder, at least as far as elec- 
tron transport was concerned). It took almost two decades, using renormalization 
group (RG) ideas, to establish that in two dimensions^ all states were again local- 
ized with arbitrary small disorder in the case of pure potential scattering, like the 
one considered by Anderson in his 1958 paper. In hindsight, given the situation in 
lower dimensions, it seems lucky that we live in a three-dimensional world, where 
many of the properties of crystalline materials we love and use are robust, at least 
to small amounts of imperfections. 

Following Anderson's stylized model, Lifshitz^ considered a somewhat more re- 
alistic model for problem of impurity bands in doped semiconductors, in which 
positional randomness of the dopants was the dominant source of disorder. Here 
the dopant atom, often less than one part in a million, dominates all the character- 
istics of the material that give it the unique properties that led to the semiconductor 
revolution. It is ironic to note that while solid state physicists studied doped semi- 
conductors quite extensively since the 1950s, nobody seems to have thought of that 
(i.e., physics of dopants) as a rare fluctuation effect! 

In modern parlance, Lifshitz considered the problem with purely off-diagonal 
(hopping) disorder, whereas Anderson had dealt with purely diagonal (onsite en- 
ergy) disorder. Many short and intermediate distance properties are better captured 
by Lifshitz's model of the impurity band, especially at low dopant densities. In his 
analysis, Lifshitz found that band tails were formed as a result of resonance between 
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pairs of sites that happened to be much closer than average, a rare fluctuation effect. 
An analogous study of band tails in the Anderson model can also be performed, 
leading to Lifshitz tails on rare configurations of clusters consisting of resonant sites, 
with many analytic results being arrived at near the band edge^ - ^. Soon thereafter, 
band tails were studied in the high-density limit by Halperin and Lax using mini- 
mum counting methods^. 

The next major appearance of rare fluctuation effects in condensed matter ap- 
pears in two studies of many-body spin systems, by McCoy and collaborators^"^ 
and by Griffitha^-l These studies, which appeared around 1968-69, dealt with clas- 
sical, Ising models of magnetism with disorder. The former was a study of a two- 
dimensional square Ising ferromagnet with uniform vertical bonds and random 
bonds in the horizontal direction, which were perfectly correlated in the vertical 
dimension - it thus appeared as a contrived model in the classical context. However, 
the results obtained were quite bizarre - it was found that the magnetic suscepti- 
bility diverged due to rare configurations of strong bonds in a finite region in the 
paramagnetic phase above the true thermodynamic phase transition temperature. 
Griffiths, on the other hand, was studying an Ising ferromagnet in any dimension 
with randomly diluted bonds. He showed that the existence of rare configuration 
of clusters with fewer than average missing bonds over finite length scales led to 
essential singularities in the thermodynamics of the system in the paramagnetic 
phase, between the transition temperature of the pure, undiluted system, and the 
actual transition temperature of the diluted system. While these weak singularities 
in the thermodynamics were experimentally undetectable, it was shown several years 
lateral that there was a concomitant slow, non-exponential decay in the long time 
dynamics of such systems. Such effects were also claimed^ to affect the long time 
dynamics of Ising spin glasses (where the bonds are ferromagnetic and antiferro- 
magnetic at random with equal probability); while non-exponential relaxatiorP^HSl 
has been a hallmark of spin glasses, whether the long time behavior is dominated 
by rare configurations remains an open issue at present. 

Roughly a decade after Griffiths and McCoy, the problem of randomly positioned 
donors surfaced agairf^, this time in the problem of the magnetic ground state of 
n-doped semiconductors, which could be modeled as a quantum spin-1/2 Heisen- 
berg antiferromagnet with an extremely high degree of disorder (large randomness) . 
This was addressed using an RG scheme that made use of large disorder, leading 
to bond distributions that were wide on a logarithmic scale 1 -. A one-dimensional 
analog of the same problem allowed an analytic treatment of the thermodynam- 
ics to be obtained in the low temperature limitP^, which was later proved to be 
asymptotically exaclPR The RG scheme showed how, in the limit of large disorder, 
the quantum system's behavior exhibits special characteristics that come about es- 
sentially because of enslavement of weak couplings by the strong ones, to a degree 
that is just not present in classical systems^. Such strong effects, which were fur- 
ther amenable to analytic RG treatments, rekindled the interest in rare fluctuation 
effects in quantum models with disorder. 
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It was recognized that the McCoy- Wu model was quite natural in the context of 
the one-dimensional random quantum Ising ferromagnet in a transverse magnetic 
fielcP3, the quantum mechanics entering through the non-commuting spin opera- 
tors (Sz and Sx) along the Ising coupling (z) and the field (x) directions. The path 
integral representation of the quantum model led to a classical model in one higher 
dimension (i.e., two), and perfect correlations between the random bonds in one 
of the dimensions represented simply the time evolution in the path integral rep- 
resentation. It was shown that approaching the ordered phase at zero temperature 
at low magnetic fields from the paramagnetic phase at high magnetic fields, the 
magnetic susceptibility diverged before the (quantum) phase transition, because of 
the effect of rare configurations of strong bonds. This generated a flurry of activity 
on transverse field quantum Ising model J^^, including those with random bonds 
(e.g. the Ising spin glass in a transverse field). At the same time, it was recognized 
that similar effects were present in random, dimerized, quantum antiferromagnetic 
spin chains^. 

The quantum spin glass models in d — 2 and higher were not amenable to 
analytic study, but Monte Carlo simulations^ showed that the divergences of the 
non-linear (spin glass) susceptibility in the paramagnetic phase due to rare config- 
urations continued to occur in d = 2 and d = 3 spatial dimensions. This strong 
effect in the thermodynamics of the quantum system, as opposed to rather weak 
effects in the classical system, can be physically understood by considering the path 
integral mapping of the quantum into a (d+1) dimensional classical system: while 
the classical random bond model has defects that are point-like in character, the 
quantum version has line defects (because of the perfect bond-correlation along the 
time direction). 

While these studies showed quite convincingly the overall significant effect of 
rare configurations for the quantum Ising spin glass, it was after a very detailed 
numerical effort involving large amounts of computer time. Looking at details as a 
function of size or geometry of clusters which taken together give rise to the singular 
behavior was not deemed possible; this situation remains true today. 

In the intervening years, there was not much further investigation of Lifshitz- 
taiP - ^ ideas in the context of the Anderson modeP; rather the excitement was 
concentrated mainly on the original issue raised by Anderson, namely, the extended- 
localized state (or the metal-insulator) transition. Starting with the pioneering nu- 
merical efforts of Kramer and MacKinnon^ and of McMillan and coworkers^, the 
study of the metal-insulator transition for non-interacting electrons in a random 
potential with different symmetry properties (orthogonal, as in potential scatter- 
ing, unitary as with a magnetic field; or symplectic as with spin-orbit scattering), 
became a virtual industry, with the most complete work done by Slevin, Ohtsuki 
and coworkers I22HH1 g ft was a complete surprise when our recent investigational 
pointed out some rather remarkable, singular behavior in Anderson's original model, 
which had escaped attention for fifty years since it had been proposed. We elaborate 
upon our findings in the next section. 
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2. Anderson's 1958 Model Revisited 

Anderson^ considered the one-electron Hamiltonian: 

i <ij> 

where i and j denote sites on a simple cubic lattice, are independent random 
variables drawn from a uniform probability distribution 

P u (e) = 1/W for W/2 < e < W/2, and zero otherwise, (2) 

and sites i and j within the angular brackets < ij > are nearest neighbors. We will 
generalize our study slightly, by considering some other probability distributions 
which are symmetric around zero (see below) and including dimensions d = 1, 2 
and 3, but restricting our attention to d-dimensional hypercubic lattices. The model 
is characterized by the dimensionless disorder, W/V. 

In three dimensions, for W = 0, the density of states (DOS) is the familiar DOS 
of the nearest neighbor tight-binding model on a simple cubic lattice extending from 
E = —6V to +6V; needless to say, all electronic states are extended, being of the 
Bloch form. As W/V is increased from zero, states near the band edges get localized, 
while those near the band center remain extended. Increasing W/V further pushes 
the mobility edge [at ±E C (W)} separating the localized tail from the extended center 
closer to the band center at E — 0, until at W ~ 16. 5V, i.e. for a disorder width of 
order the original bandwidth (12V"), all states become localized [6]. In d = 1 and 2, 
as discussed earlier, all states become localized for any nonzero W/V. 

The other distributions we will discuss are Gaussian: 

P G (e) = {2^ 2 )- 1 / 2 exp{-e 2 /2a 2 ), (3) 

where a — W/\/V2 gives the same variance as Pjj, the triangular distribution: 

P T (e) = 2(1 - 2\e\/W)/W for |e| < W/2, ), otherwise, (4) 

as well as a distribution with a power-law behavior at the edges (±W/2) and zero 
density of states at the origin: 

P PL (e) cx \e\/(W 2 - 4e 2 ) Q ~ 1 for |e| < W/2, otherwise (5) 

In Eq. [5] a > 0, so that the distribution is integrable. 

Fig. [T] depicts the first three of these distributions of the on-site energy, given 
by Eq. [2]|4j These distributions will help illustrate our findings and the conditions 
under which we get them. 

In numerical studies of Anderson localization, two main quantities that have 
proved useful to calculate the localization length at any energy E, £(E), have been 
the inverse participation rations (IPR(E)) , and the Lyapunov exponent!^ Ly(E) . 
The former is easily defined for any wavef unction |\& >= J2i c iV >! as 



/pp* = (]T c f)/(E c ') 2 - 

i i 



(6) 
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Then 

(IPR(E)) = (IPR 9 6(E - En,)), (7) 

where the average is over states \& that are at an energy in a small window around E, 
which can be made arbitrarily small in the thermodynamic limit, when the typical 
energy spacing scales as 1/L d , where L is the linear dimension of the system. For a 
state that is localized on a length scale £, it is easily seen that (I PR) ~ C d - F° r an 
extended state, on the other hand, (IPR) ~ L~ d , where L is the linear dimension 
of the system. 

The Lyapunov exponent at energy E is computed for systems of finite cross- 
section (M d ~ 1 ) and long length in one dimension, which is taken to the thermody- 
namic limit. It is defined in terms of the eigenvalues of the transfer matrix connecting 
the amplitudes of a wavefunction at position i and (i + 1) with those at (i + 1) and 
(i + 2). In dimensions d > 1, one has to evaluate it for different M and take the limit 
as M oo. For our purpose, it suffices to note that the limit, Ly(E) characterizes 
the exponential decay of the wavefunction at long distances, at an energy E, and 
therefore, 



(8) 

[Ly(E)] d , 



Ly(E) ~ l/£(E) 

On the basis of ^ and one expects therefore that (IPR(E)) 
and in particular, for d = 1, they should be proportional to each other. 

In the next section, we describe the results of our numerical investigations, by 
using d = 1, where states are localized at all energies for nonzero disorder (W/V > 
0). We discuss the situation in higher dimensions briefly in Section 4. 



3. Numerical Results (d=l) 

Fig. [2] shows (IPR(E)) and Ly(E) as a function of E for the Anderson Model (Eq. 
1) with a uniform distribution (Eq. 2), in d = 1 for W/V = 6. Since all states 
are localized, results converge rapidly as a function of size, and do not require 
extensive numerical effort. As can be seen, starting from the band center, both 
(IPR(E)) and Ly(E) increase with \E\, indicating a decrease in the localization 
length with increasing \E\, as may be expected (band tails localize first even in 
d = 3). However, as one goes to larger \E\, Ly(E) continues to rise, but (IPR(E)) 
goes through a maximum, and decreases, apparently going to zero at the band 

P{a) 




-W/2 W/2 -w/2 w/2 

Fig. 1. From left to right, uniform, Gaussian and triangular distribution of disorder. 
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edge. This is consistent with Lifshitz's considerations — the states at the absolute 
band edge for finite W [given by Eg = ±(W/2 + 2dV) for hypercubic lattices 
in d-dimensions] are due to states on rare clusters of sites with on-site energies 
all at the edge of the disorder probability distribution, i.e. for ej = ±W/2 for all 
sites. Of course, such a configuration is exponentially rare in the size of the cluster 
~ exp[— cL d ] where L is the linear dimension of the cluster, and has a vanishing 
probability in the thermodynamic limit. Of course, in such a case, the eigenstates 
are spread all over the cluster since there is no disorder, and would therefore have an 
I PR ~ 1/L d . Particle in a box considerations suggest that the energy of these states 
would be ~ 1/L 2 away from the true band edge, implying that in the asymptotic 
region, (IPR(E)) ~ \E B ~ E\ d /' 2 , which -> as E -t E B . Clearly in this (Lifshitz) 
regime, the extent of the wavefunction (measured by IPR) and the exponential decay 
length [] measured by Ly(E)] are not the same; one has two length scales describing 
the electronic wavefunction. 

What is more striking than even the downturn is the manner in which the IPR 
changes course: there is a clear cusp in the data, implying non-analytic behavior. 
Furthermore, this non-analytic downturn happens at an energy Er, where the den- 
sity of states (DOS) is a maximum, not deep in the exponential tail of the DOS, 
as shown in Fig. [3] In fact, the DOS also displays a cusp at the same energy. (This 
does not contradict the result of Edwards and Thouless^ who showed that the 
DOS is analytic in an energy interval (W/2 — 2dV) around the band center; the 
singular behavior we see is outside this region). The cusp- like behavior is seen for 
intermediate and large (W/V); in d = 1, we see it clearly for all W/V > 3.8. We 
remark that Ly(E) does not show any sign of the singular behaviour at Er. 




Energy 

Fig. 2. (IPR(E)) & Ly(E) for the uniform distribution with W/V = 6. 
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Fig. 3. Density of states for the uniform distribution with W/V = 6. 
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As explained in Ref. this singularity appears to be a consequence of the 
fact that typical Anderson localized states with a central site do not exist beyond 
Er] the only states that are present beyond it are states that involve resonant 
tunneling between two or more nearest neighbor sites. In the Anderson model, for 
large disorder W, typical Anderson localized states have their energies perturbed 
by 0(V 2 /W), whereas resonant states on two or more sites have energies shifted 
by 0(V); increasing the number of resonant sites only changes the prefactor of 
the 0(V) term, not the functional dependence. As a result, resonant states on 
sites near the disorder band edge pull out of the quasi-continuum of the typical 
Anderson-localized states, leaving behind non-analyticity in the density and nature 
of electronic states at the edge. 

Er thus appears to be the beginning of the Lifshitz tail, whose precise nature 
depends on short-distance rather than long-distance physics. Furthermore, since 
the smallest resonant cluster consists of 2 sites, its probability is not exponentially 
small, but only 0(V/W) for large W, so the number of resonant states beyond Er 
can be significant. For the case of the uniform distribution in d = 1, we finoP that 
the maximum percentage of states beyond Er is « 17% (for W/V rj 3.8). This 
fraction of tail states can hardly be called rare; hence the quotation marks in the 
title of this paper. We note in passing that we have shown [30] by explicit analytic 
calculation for a two-site Anderson model, that such cusp-like singularities exist in 
both (IPR(E)) and DOS(E). 

The non-analytic behavior characterizing the separation of resonant states from 
typical Anderson localized states appears to exist for all forms of bounded disorder. 
Fig. 4 shows (I PR) and DOS as a function of E for the triangular distribution 
(Eq. 4) for W/V = 6. Again the sharp maximum in (I PR) occurs at a DOS that is 
not exponentially small, despite the smoother edge to the probability distribution 
of the disorder. On the other hand, for the Gaussian distribution of disorder (Eq. 
[3]), we find that both Ly(E) and (IPR(E)) increase monotonically with \E\, with 
no sign of a turnaround or strong singularity for up to values of \E\ where the DOS 
is down to 10~ 4 of its value at the band center. This is consistent with Wegner's 
result!^ that the DOS is analytic at all E for Gaussian disorder. 

For distributions like that in Eq. [5] with several distinct edges (\E\ = W/2, and 
E = 0), we find that the system develops more singular points dividing typical and 




Fig. 4. DOS and (IPR(E)) for the triangular distribution of onsite energies with W/V = 6. 
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resonant eigenstates. This is currently under further investigation. 

On the basis of the above findings, we conclude that the Anderson model with 
bounded disorder in d — 1, when all states are localized, has a critical energy En 
at least at large disorder, which demarcates the transition from typical Anderson 
localized states with a central site, to states that necessarily involve resonance be- 
tween two or more nearest neighbor sites. This is depicted in Fig. [5j When these 
resonant Lifshitz-like states occur at an energy that is forbidden to typical Ander- 
son localized states, they maintain their strongly resonant character. On the other 
hand, when they occur at energies in the range allowed to typical Anderson localized 
states, they mix with the more abundant typical states, and lose their special reso- 
nant character. In that sense, this phenomenon appears to be analogous to the fact 
that localized states can only exist at energies where extended states do not; where 
they do, they mix and lose their localized character. Whether there is a further 
hierarchy among resonant states on different number of sites, and whether these 
are accompanied by (presumably weaker) singularities or not, remains a topic for 
further investigation. 

4. Higher Dimensions and Concluding Remarks 

In conclusion, we have shown that Anderson's original model of localization exhibits 
behavior that is richer than believed, so far. for the uniform distribution, as well as 
other bounded distributions of diagonal disorder, in the localized phase, the inverse 
participating ratio characterizing the extent of localized states exhibits a distinct, 
non-analytic behavior, which is to some extent reflected in the density of states as 
well. The energy En at which this non-analyticity occurs separates typical Anderson 
localized states from states that have a resonant character on various length scales. 

This phenomenon, described above in Section 3 for d = 1 is not special to one 
dimension. As shown in Ref. 1291 cusp-like singularities exist for large disorder in 
d = 2, and also in d = 3, in the insulating phase. Of course, detailed numbers are 
dimension dependent. It is interesting to speculate whether in d = 3, the existence 
of a second critical energy Eq denoting the transition from extended to localized 
behavior will lead to an interplay of some sort. For example, would \Er\ always 
be bigger than \Eq\ for all W, or could there be a transition from extended states 
directly into localized, resonant states as a function of energy for some value of the 



N(E) 



Resonant 
Lifshitz — like 

states 




Resonant 
Lifshitz — like 

states 



Fig. 5. A schematic illustration showing the occurence of different kinds of wavefunctions in the 
insulating phase. 



March 2, 2013 



10 R.N. Bhatt, Sonika Johri 

disorder? If so, would that be any different from the standard localization transition? 
In this regard, it is interesting to note that numerical investigations^^!]!!!] of the 
localization transition in the Anderson model in d = 3 have found universal behavior 
for the transition at the band center [E = 0) as a function of disorder (W), but 
have found non-universal results with similar size systems, when they study the 
transition at fixed W as a function of energy (E). It is tempting to speculate that 
this may be connected to the change in the nature of the localized state as a function 
of energy at -Br, which is quite distinct from the true band edge Eb- Clearly, further 
research needs to be done to determine if the speculation has merit. 

Another interesting result of our findings is that the original Anderson modeP 
provides a platform for an in-depth study of rare fluctuation physics. While ana- 
lytic methods are most useful in the asymptotic exponential tail, where numerical 
methods are of little use, numerical methods can help provide detailed results due 
to resonant clusters of small finite sizes. If analytic methods can help decipher the 
nature of further (weaker) singularities due to the transitions from n-site clusters to 
(n + l)-site clusters that may exist as one descends down the tail of the distribution, 
there is the possibility that numerical methods may be able to pick them up quan- 
titatively. This is because Anderson localization, being a single particle problem, 
requires time that grows as a power of the system size. Thus it is amenable to a 
much more intensive numerical investigation than many-body Hamiltonians with 
exponential growth of the Hilbert space with size. In addition, since most of the 
physics lies in the insulating phase, the sizes that will need to be considered will 
likely also be manageable. 
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